(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(11) EP 0 664 037 B1 

EUROPEAN PATENT SPECIFICATION 



(45) Date of publication and mention 
of the grant of the patent: 
19.12.2001 Bulletin 2001/51 

(21) Application number: 92921735.4 

(22) Date of filing: 10.10.1992 



(51) IntClJ: G06K 9/00 

(86) International application number: 
PCT/US92/08673 

(87) International publication number: 

WO 94/09446 (28.04.1994 Gazette 1994/10) 



(54) BIOMETRIC PERSONAL IDENTIFICATION SYSTEM BASED ON IRIS ANALYSIS 

BIOMETRISCHES PERSONENIDENTIFIZIERUNGSSYSTEM AUF DER BASIS VON 
IRIS-ANALYSE 

SYSTEME D'IDENTIFICATION DE PERSONNE BIOMETRIQUE BASE SUR L'ANALYSE DE L'IRIS 



(84) Designated Contracting States: 

AT BE CH DE DK ES FR GB GR IE IT LI LU MC NL 
SE 

(43) Date of publication of application: 
26.07.1995 Bulletin 1995/30 

(73) Proprietor: Iridian Technologies, Inc. 
Moorestown, New Jersey 08057 (US) 

(72) Inventor: Daugman, John G. 

Huntingdon, Cambridgeshire PE18 9JX (GB) 

(74) Representative: 

Robinson, Nigel Alexander Julian 
D. Young & Co., 21 New Fetter Lane 
London EC4A 1DA (GB) 



QQ 

CO 

o 

CO 
CD 



Q. 
LU 



(56) References cited: 
US-A-4 620 318 
US-A-5 016 282 



US-A- 4 641 349 



• PROCEEDINGS. THE INSTITUTE OF 

ELECTRICAL AND ELECTRONICS ENGINEERS 
1992 INTERNATIONAL CARNAHAN 
CONFERENCE ON SECURITY TECHNOLOGY: 
CRIME COUNTERMEASURES (CAT. 
NO.CH31 19-5/92), ATLANTA, GA, USA, 14-16 
OCT. 1992, ISBN 0-7803-0568-X, 1992, NEW 
YORK, NY, USA, IEEE, USA, pages 50-60, 
DAUGMAN J 'High confidence personal 
identification by rapid video analysis of iris 
texture 1 



Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give 
notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in 
a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 
99(1) European Patent Convention). 



Printed by Jouve, 75001 PARIS (FR) 



EP 0 664 037 B1 



Description 

[0001] This invention is related to the field of personal identification, and more specifically to the field of automated 
identification of humans by biometric indicia. 

5 [0002] Identification of humans is a goal as ancient as humanity itself. As technology and services have developed 
in the modern world, human activities and transactions have proliferated in which rapid and reliable personal identifi- 
cation is required. Examples include passport control, computer login control, bank automatic teller machines and other 
transactions authorization, premises access control, and security systems generally. All such identification efforts share 
the common goals of speed, reliability, and automation. 

w [0003] The use of biometric indicia for identification purposes requires that a particular biometric factor be unique 
for each individual, that it be readily measured, and that it be invariant over time. Although many indicia have been 
proposed over the years, fingerprints are perhaps the most familiar example of a successful biometric identification 
scheme. As is well known, no two fingerprints are the same, and they do not change except through injury or surgery, 
it is equally clear, however, that identification through fingerprints suffers from the significant drawback of requiring 

15 physical contact with the person. No method exists for obtaining a fingerprint from a distance, not does any such method 
appear likely. 

[0004] A biometric indicator that has been largely ignored by the art is the iris. The iris of every human eye has a 
unique texture of high complexity, which proves to be essentially immutable over a person's life. No two irises are 
identical in texture or detail, even in the same person. As an internal organ of the eye the iris is well protected from the 

20 external environment, yet it is easily visible even from yards away as a colored disk, behind the clear protective window 
of the eye's cornea, surrounded by the white tissue of the eye. Although the iris stretches and contracts to adjust the 
size of the pupil in response to light, its detailed texture remains largely unaltered apart from stretching and shrinking. 
Such distortions in the texture can readily be reversed mathematically in analyzing an iris image, to extract and encode 
an iris signature that remains the same over a wide range of pupillary dilations. The richness, uniqueness, and immu- 

25 tability of iris texture, as well as its external visibility, make the iris suitable for automated and highly reliable personal 
identification. The registration and identification of the iris can be performed using a video camera without any physical 
contact, automatically and unobtrusively. 

[0005] By comparison, other biometrics such as signatures, photographs, fingerprints, voice prints, and retinal blood 
vessel patterns all have significant drawbacks. Although signatures and photographs are cheap and easy to obtain 
30 and store, they are impossible to identify automatically with assurance, and are easily forged. Electronically recorded 
voice prints are susceptible to changes in a person's voice, and they can be counterfeited. Fingerprints or hand prints 
require physical contact, and they also can be counterfeited and marred by artifacts. 

[0006] Iris identification is not to be confused with retinal identification. The iris is easy to see and can readily be 
imaged with a video camera. The retina, on the contrary, is hidden deep within the eye, and is difficult to see. Common 
35 conditions such as small pupils or cataracts make it difficult or impossible to see the retina, but they do not affect the 
visibility of the iris. 

[0007] The only previous attempt to take advantage of these favorable characteristics of the iris for a personal iden- 
tification system is seen in U.S. Patent No. 4,641,349, issued to Flom and Safir and entitled "Iris Recognition System, 
" (hereinafter '349). The '349 reference discloses the general concept of using the iris as a method of identification, 

*o but it does not describe a developed embodiment of such a system. It does not disclose automatic means to find and 
isolate the iris within an image, regardless of the location and size of the iris, nor means to extract and encode its 
texture. The '349 reference uses adaptive illumination to bring the pupil to a predetermined size in order to overcome 
the problems of variation in iris pattern with variation in pupil size. Moreover, the '349 reference does not describe any 
method for computing an identification decision once a list of features has been compiled. In general, a listing of features 

45 from two iris images will partially agree and partially disagree, whether or not they originated from the same iris. Gen- 
erally such lists will also differ in the number of features they comprise. No theoretical or mathematical formulation was 
provided for basing decisions on such comparisons between incommensurate data sets. Moreover, no method was 
disclosed for calculating the confidence levels associated with identifications. 

[0008] The broad object of the present invention is to provide a method and apparatus for identifying persons based 
50 on the iris of either eye. 

[0009] A further object of the invention is to allow extremely reliable and rapid identification of a person, in approxi- 
mately one second. 

[001 0] Another object of preferred embodiments of the invention is to provide an identification system that calculates 
confidence levels for any identification decision, on an objective and rigorous basis. 
55 [001 1] A still further object of the invention is to provide an identification system that provides identification without 
action from the subject and without making physical contact with the subject. 

[0012] These and other objects are achieved in the present invention by a method and apparatus for uniquely iden- 
tifying a particular human being as defined in claims 1 and 24. The following steps are carried out in detail: first, the 



2 



EP 0 664 037 B1 



system acquires through a video camera a digitized image of an eye of the human to be identified. Then, it isolates 
the iris if it is present within the image and defines a circular pupillary boundary between the iris and pupil portions of 
the image, and it defines another circular boundary between the iris and sclera portions of the image, using arcs that 
are not necessarily concentric with the pupillary boundary. The system of the invention establishes a polar coordinate 

5 system on the isolated iris image, the origin of the coordinate system being the center of the circular pupillary boundary. 
It then defines a plurality of annular analysis bands within the iris image, these analysis bands excluding certain prese- 
lected portions of the iris image likely to be occluded by the eyelids, eyelashes, or specular reflection from an illuminator. 
[0013] The portion of the iris image lying within these annular analysis bands is analyzed and encoded employing a 
special signal processing means comprising a multi-scale, self-similar set of quadrature bandpass filters in polar co- 

10 ordinates, to generate an iris code of fixed length and having a universal format for all irises. The resulting code is 
stored as a reference code. Because of the universal format and length of all such iris codes, comparisons among 
different iris codes are extremely efficient and simple. Specifically, a comparison between any two iris codes is achieved 
by computing the elementary logical XOR (exclusive-OR logical operation) between all their corresponding bits, and 
then computing the norm of the resulting binary vector. This comparison measure can also be described as the Ham- 

15 ming distance between the two iris code vectors. The universal format of iris codes also lends itself to rapid parallel 
search across large data bases of stored reference iris codes in order to determine the identity of an individual. 
[0014] To perform an identification using the reference code, the system generates from an identification subject an 
identification code according to the providing, defining, determining, establishing, describing and analyzing steps. Then, 
the system compares the identification code with the reference code, to ascertain the Hamming distance between the 

20 codes. This distance is then converted into a calculated likelihood that the two codes originated from the same iris, 
and hence from the same person, by computing the probability that the observed matching fraction of bits in the two 
codes could match by chance if the two codes were independent. A preselected criterion applied to this measured 
Hamming distance generates a "yes" or "no" decision, and the confidence level for the decision is provided by the 
calculated probability. 

25 [001 5] These together with other objects of the invention, along with various features of novelty which characterize 
the invention, are pointed out with particularity in the claims annexed hereto and forming a part of this disclosure. For 
a better understanding of the invention, its operating advantages and the specific objects attained by its uses, reference 
should be had to the accompanying drawings and descriptive matter in which there is illustrated a preferred embodiment 
of the invention. 

30 

Fig. 1 is a block diagram depicting the principal steps of the process of the present invention; 

Fig. 2 is a photograph of a human eye, overlaid to demonstrate the process of the present invention; 

35 Fig. 3 illustrates the quadrature bandpass filters employed by the present invention as image convolution kernels 

to extract iris structure at many scales of analysis. 

Figs. 4A-4C are a set of plots illustrating how iris image data is converted into iris code bits by the quadrature 
bandpass filters according to the present invention; 

40 

Fig. 5 is a plot showing the Hamming distances for "imposters", i.e. , for comparisons between iris codes computed 
from different irises; 

Fig. 6 is a plot showing formulation for statistical decision theory; 

45 

Fig. 7 is a plot illustrating the theoretical effects of manipulation of decision criterion; 

Fig. 8 is a plot showing the Hamming distances for "authentics", i.e., for comparisons between iris codes computed 
on different occasions from the same iris; 

50 

Fig. 9 is a plot showing Hamming distances for authentics and imposters, combined; 
Fig. 10 is a plot showing binomial fit to Hamming distance distributions; 
55 Fig. 11 is a plot of error as a function of Hamming distance criterion; and 

Fig. 12 is a table of performance rates achieved with the present invention. 
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[0016] An embodiment of the present invention is shown in schematic form in Figure 1 and comprises a block diagram 
depicting the principal steps in developing an iris identification code for a person, and then using that code to make 
an identification. The process will be discussed in overall terms, followed by a detailed analysis. 
[0017] The iris of the human eye is a complex structure comprising muscle, a trabecular meshwork of pectinate 

5 ligament, connective tissue, blood vessels, and chromatophores. Externally it presents a visible texture with both radial 
and angular variation arising from contraction furrows, collagenous stromal fibers, filaments, serpentine vasculature, 
rings, crypts, and freckles; taken altogether, these constitute a distinctive "fingerprint." The magnified optical image of 
a human iris thus constitutes a plausible biometric signature for establishing or confirming personal identity. Further 
properties of the iris that lend themselves to this purpose, and render it potentially superior to fingerprints for automatic 

10 identification systems, include the impossibility of surgically modifying its texture without unacceptable risk; its inherent 
protection and isolation from the physical environment, being an internal organ of the eye, behind the cornea and 
aqueous humour; and its easily monitored physiological response to light, which provides a natural test against artifice. 
Additional technical advantages over fingerprints for automatic recognition systems include the ease of registering the 
iris optically without physical contact, and the intrinsic polar geometry of the iris, which imparts a natural coordinate 

15 system and origin. 

[0018] Unknown until the research leading to the present invention was whether there are sufficient degrees-of- 
freedom, or variation in the iris across individuals, to impart to it the same singularity as a conventional fingerprint. Also 
uncertain was whether efficient algorithms could be developed to extract detailed iris structure reliably from a video 
image, generate a compact code (of minuscule length compared with image data size), and render a decision about 
20 identity with high statistical confidence, all within less than one second of processing time on conventional equipment. 
The present invention resolves all of these questions affirmatively. 

[0019] At the broadest level, the system of the present invention can be envisioned in five steps. First, an image of 
the eye to be analyzed must be acquired in digital form suitable for analysis, as shown in block 10 of Figure 1 . Then, 
the iris portion of the image must be defined and isolated (blocks 12, 14, 16 and 18). The defined area of the image 

25 must next be analyzed to produce an iris code (block 20). It should be noted that the first iris code generated for a 
particular iris is stored as a reference code (block 22). In subsequent encounters, the system uses the reference code 
to make an identification by comparing a presented code (block 24) with the reference code to obtain a Hamming 
distance (block 26). This data permits the system to establish, confirm, or disconfirm the identity of the subject (block 
28), and to calculate a confidence level for the decision (block 30). The rapidity of the process also enables passive 

30 identification by exhaustive search through a large database, rather than merely validating a single authorization code. 
[0020] In a practical application of this system, a digitized image as illustrated in Figure 2 shows an eye 100 with an 
iris 102 surrounding a pupil 104. The eye's sclera portion 105 which is the white portion in turn surrounds the iris 102. 
The first step in processing the image is to locate the pupillary boundary 106, separating the pupil 104 from the iris 
102, to a high degree of accuracy (block 12, Fig. 1). This step is critical to insure that identical portions of the iris are 

35 assigned identical coordinates every time an image is analyzed, regardless of the degree of pupillary dilation. 

[0021] The inner boundary of the iris, forming the pupil, can be accurately determined by exploiting the fact that the 
boundary of the pupil is essentially a circular edge. As seen in Figure 2, the pupil 104 is generally dark while the iris 
102 is lighter, with varied pigmentation. However, this relationship may sometimes be reversed, for example in eyes 
with dark irises and some cloudiness of the internal lens, or because of optically co-axial illumination (directly into the 

40 eye), in which case light is reflected back from the retina and out again through the pupil. A further reason that the 
image of the pupil may be bright is because of specular reflections from the cornea. A general purpose method for 
finding the pupillary boundary should be robust enough to function reliably whether or not the region of the pupil is 
actually darker than the iris. In the present invention, a system of integrating evidence about the true pupillary boundary 
has been developed which has the desired robust behavior and accuracy. 

45 [0022] The method of the present invention detects the pupillary boundary as an abrupt and sudden change in bright- 
ness when summed along a circle whose radius is steadily increasing. This sudden change will be maximum if the 
circle has its center near the true center of the pupil, and when its radius matches the true radius of the pupil. Thus, 
the image processing problem of finding the pupil can be formulated as an optimization problem, in which a series of 
"exploding circles" (steadily increasing radii) are positioned with their center coordinates located at each one of a 

50 number of trial points on a grid. For each exploding circle, and for each value of its radius, the total image brightness 
is summed over a fixed number of points lying on this circle. (Using a constant number of points on each circle, typically 
128, avoids an automatic increase in the summed brightness simply due to increasing circumference.) The system 
searches for the maximum rate of change in this quantity as radius expands. For the candidate circle that best describes 
the pupillary boundary, there will be a sudden "spike" in the rate-of-change of luminance summed around its perimeter, 

55 when its radius just matches that of the pupillary boundary. This spike will be larger for a circle that shares the pupil's 
center coordinates and radius, than for all other circles. In this manner, the problem of precisely locating the pupil has 
been converted into an optimization problem in which a three-parameter space is searched for the best combination 
of circle center coordinate (Xq, y 0 ) and radius r. 
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[0023] The process of the instant invention can be described mathematically as summing a contour integral of image 
intensity, l(x,y), over the arc (ds) of a circle whose radius is r and whose center coordinates are (x Q , y Q ); and then 
computing the partial derivative of this quantity with respect to r as the radius increases. The maximum absolute value 
of this derivative is sought, over the space of the three parameters (x 0 ,y 0 ,r): 



The partial derivative with respect to r can be a smoothed, or blurred, derivative for noise immunity, and it may also 
be converted into a percentage change (dividing by the current value of the contour integral) for enhanced noise im- 
munity. The method also has some intrinsic noise immunity, because the contour integral inherently integrates data 

15 over a contour, so any anomalous excursions in pixel brightness tend to average out. 

[0024] The search process in the three-parameter space is directed by gradient ascent, or "hill-climbing." If a can- 
didate series of exploding circles is partially within the pupil, the value of the quantity defined in Equation (1) will be 
larger than for other circles that are not. The closer the center of the concentric circles comes to the true pupillary 
center, the larger will be this quantity. Similarly, the quantity within Equation (1) will be larger for circles of appropriate 

20 radius. Thus, the method can find the optimal combination of the three parameters by an iterative search process, in 
which the step size of changes in the three parameters is decreasing with each successive iteration. By moving always 
in the best direction (in which the rate of improvement is greatest) in the three-parameter space, and by taking stepsizes 
proportional to the rate of improvement, and by decreasing these stepsizes on each successive iteration, the method 
rapidly converges. After usually only four or five iterations, the optimal values of the three parameters have been 

25 determined to within less than one pixel. The values for (x 0 ,y Q , r) determine the estimated pupillary boundary, as well 
as the origin of the polar coordinate system for subsequent iris analysis. 

[0025] This efficient method of finding and tracking the pupillary boundary also provides an important safeguard 
against imposters. One obvious method for trying to defeat an identification system based on iris patterns would be to 
present to the video camera a photograph of anothers person's eye, or even to wear contact lenses imprinted with the 

30 image of an authorized iris. However, an important feature of a living eye is that the pupil diameter continuously un- 
dergoes small oscillations ("hippus") once or twice per second, even under uniform lighting. A photograph of an iris, 
or a contact lens imprinted with an iris image, would not exhibit such variation in time. Because the process described 
above for finding and tracking the pupillary boundary is so rapid, it is possible to acquire several images in succession 
and to monitor the pupillary diameter over time. An absence of hippus oscillations or other small variations in iris pattern 

35 over time would constitute evidence that a photograph or simulacrum were being presented, rather than a living iris, 
and this would indicate an intrusion attempt. This ability to discriminate between a living iris and a simulacrum or 
photograph is an important security asset, made possible by the rapid means for defining and tracking the pupillary 
boundary. 

[0026] Once the boundary and center have been determined, the next step is to locate the outer boundary of the 

40 iris, or limbus, where it meets the sclera. An important consideration here is the fact that the pupil is not always centered 
within the iris. The radial distances to the right and left limbus may vary by as much as 20%, and hence both distances 
must be computed in order to generate an appropriate iris coordinate system. A further consideration is the fact that 
the upper and lower eyelids generally obscure the top and bottom boundaries of the iris, and hence these regions 
should be excluded from iris analysis. 

45 [0027] The same general method of "exploding circles" that yielded accurate determination of the pupillary boundary 
can be used for finding the outer boundary of the iris, but with two modifications. First, given the upper and lower eyelid 
occlusions and the generally unequal left and right limbic distances, the method is restricted to just two arcs along the 
horizontal meridian, one at 3-o'clock and one at 9-o'clock, each subtending 7t/4 radians (45 degrees). The distances 
to these two boundaries on either side of the iris are measured separately. Second, because of possible concentric 

50 texture within the iris that could produce a maximum in Equation (1), the contour integral used earlier for the pupillary 
boundary is replaced by an area integral that blurs outs iris detail in searching for the limbus. In effect, the "exploding 
circles" are replaced by two horizontally "exploding pie wedges," that search for a sustained luminance step signifying 
the sclera on either side. As before, the search process remains one of finding the maximum in the rate-of-change of 
integrated luminance as the radius of expansion increases. After compensation for the increasing area of luminance 

55 integration, the maximum in this derivative with respect to radius invariably corresponds to the correct left and right 
boundaries of the iris. 

[0028] Mathematically, this operation is implemented by searching for the value of r (the distance from the pupillary 
center to either the right or the left) which maximizes the expression 
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max rell . 5roil0ro) - j p=r 6 — J fm4 _ 9/t J(p, e)pdpdO 



where r Q is the pupillary radius (computed earlier), 5 is a thin radial shell (typically 0.1 r 0 ), I (p,8) is the image intensity, 
now expressed in terms of polar coordinates p and 6, and <(> equals either 0 or n t corresponding to either the 3-o'clock 

10 or the 9-o'clock meridian, respectively. It has been found successful to compute this expression for values of r between 
1.5r c , and 10r o (i.e., 1.5 times to 10 times the pupillary radius) in searching for the outer boundaries of the iris, thus 
covering a wide range of possible relative diameters of pupil and iris. Similarly, the choice of +/- ji/8 radians as the arc 
angle of integration in Equation (2) has proven a useful horizontal angular delimiter for the pie wedges, to avoid the 
upper and lower eyelids. The calculated results from Equation (2) are shown in Figure 2 as the series of white dots, 

75 110-1 and 110-r on the iris 102, which correctly correspond to the left and right boundaries of the iris. 

[0029] In summary Equation (1) finds the inner boundary of the iris, i.e., the pupillary boundary. This equation gen- 
erates a series of "exploding circles" at various center positions, searching iteratively for the one combination of pa- 
rameters (x Q and y Q center, and radius r) at which the integrated luminance along the circle undergoes the greatest 
absolute value of rate-of-change. Hence, we seek the maximum absolute value of the partial derivative with respect 

20 to r, of the contour integral of luminance along the circle. This search covers the x 0 , y 0 , r parameter-space, in a very 
efficient iterative process of gradient-ascent. Equation (2) finds the outer boundary of the iris, namely the limbus, where 
the white sclera begins. The same process of exploding circles as in Equation (1) would work, but for (i) upper and 
lower eyelid occlusion which could cause difficulties; and (ii) the fact that the iris is less uniform than the pupil, and 
may have large "circular edges" itself that could trap the algorithm of Equation (1). So Equation (2) instead specifies 

25 a series of "exploding pie wedges" in the horizontal meridian (hence avoiding upper and lower eyelids), and integrates 
luminance within pie wedges rather than just along a circle. Therefore, Equation (2) specifies an area integral in polar 
coordinates that is differentiated with respect to radius, rather than a contour integral as in Equation (1). 
[0030] With the locations of the pupillary boundary and the limbus established, and the origin of polar coordinates 
fixed at the center of the pupil, a series of zones of analysis are assigned to regions of the iris. These are defined 

30 concentrically at fixed linear fractions of the radial distance between pupil and limbus, whatever may be the overall 
size of the iris in a given image, to achieve size invariance in the code. Thus the polar coordinate system for the iris is 
dimensionless in both its angular and radial coordinates. Since the iris can be approximately modelled as a rubber 
sheet that stretches and contracts with the pupillary reflex, its texture and markings stretch and shrink accordingly. 
These distortions are eliminated by the use of a radial coordinate that marks off distance simply as different fractions 

35 of the total distance from the inner boundary of the iris (the pupil) to its outer boundary. Thus a given iris, in different 
states of pupillary dilation on different occasions, should generate approximately the same iris code. A second purpose 
served by this dimensionless coordinate system is that differences in the overall size of the iris image per se, due to 
image acquisition from different distances, will not change the computed iris code. 

[0031] Since the pupil in general is not perfectly centered horizontally within the iris, it is necessary forthe fractionation 
40 to be based on a linear combination of right and left limbus estimates, weighted cosinusoidally by angle. Due to frequent 
partial occlusion of the upper iris by the upper eyelid, and specular reflection from the cornea obscuring part of the 
lower iris, these areas are excluded from analysis and encoding. An illustration of these ultimate zones of analysis, 
superimposed upon a particular iris image, may be seen in Figure 2. 

[0032] Specifically, the portion of the iris to be analyzed is mapped and subdivided into analysis bands 112 (see 
45 Figure 2). These analysis bands are defined in a special polar coordinate system whose radial coordinate may be 
slightly distorted if, as frequently occurs, the inner and outer boundaries of a particular iris are not concentric. Specif- 
ically, for any angular coordinate around the iris, the radial coordinate r of some point is defined by its fraction of the 
distance from the pupillary boundary to the sclera, along that ray. Thus, just as an angular coordinate is (classically) 
a dimensionless quantity between 0 degrees and 360 degrees, so also the radial coordinate is dimensionless in this 
50 bounded system, lying always in the interval between 0 and 1, regardless of the overall size of the iris image and 
regardless of the degree of pupillary dilation. Hence this doubly-dimensionless polar coordinate system is inherently 
size-invariant (thus inherently compensating for variation in the distance from the eye to the video camera). Similarly, 
the coordinate system also inherently compensates for any non-concentricity of the inner and outer boundaries of an iris. 
[0033] Four additional special features of the analysis bands 112 are needed to compensate for the departure of 
55 many iris images from the ideal, annular stereotype. First, since the pupil 104 itself often has an irregular boundary, 
the innermost analysis band starts with a radius of about 1 . 1 times the pupil's average radius in order to ensure exclusion 
of the pupil entirely. Similarly, since the transition from iris 102 to sclera 105 may likewise be irregular and non-circular, 
the outermost analysis band extends radially only to about 80% of the distance to the outer boundary of the iris (as 
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measured to the right and left, with cosinusoidal weighing in intermediate angles). Third, provision must be made for 
the occlusion of the upper and lower portions of the iris by the eyelids, and fourth, for a specular reflection that may 
cover part of the iris if an oblique source of illumination is used (typically from below). These spurious features are 
excluded by restricting the outermost analysis bands to two sectors centered on the horizontal meridian, thus avoiding 
5 the regions likely to be occluded by the upper and lower eyelids and by also excluding a narrow notch around the 
6-o'clock position for the specular reflection of illumination from below. These excluded areas are noted in Figure 2 as 
114a and 114b. It is preferred to divide the iris area into as many as eight annular bands 112 for analysis, at fixed 
radial fractions of the distance between the inner and outer boundaries as defined above. 

[0034] Having accurately defined the image area subject to analysis, the system then processes the data obtained 

10 from that area to generate the identification code, as depicted in block 20 of Figure 1 . Unlike the system described in 
the prior art, the present invention does not depend on controlling the amount of pupillary dilation. Rather, because of 
the dimensionless radial coordinate that simply measures off certain fractions of the distance from the inner boundary 
to the outer boundary of the iris, any given piece of iris tissue will always fall into the same zone of analysis, with the 
same position coordinates, regardless of how stretched or unstretched the iris happens to be due to pupillary dilation. 

15 This dimensionless coordinate system exploits the fact that the stretching of the iris can be approximated as the stretch- 
ing of a rubber sheet, so that its markings can still be recovered mathematically in undistorted form because the co- 
ordinate system is stretched an equal amount. Therefore, the texture of the iris is always encoded into essentially the 
same iris code, regardless of the degree of pupillary dilation, and regardless also of the overall size of the iris image. 
[0035] An effective strategy for extracting textural information from images, such as the detailed patterns of the iris, 

20 is convolution with quadrature bandpass filters such as 2-D Gabor filters. These 2-D filters were proposed by the 
inventor in 1980 and 1985 both as a framework for understanding orientation- and frequency-selective receptive field 
properties observed in neurons in the primate visual cortex, and as useful operators for practical image analysis prob- 
lems. See Daugman, J. (1980) "Two-Dimensional Spectral Analysis Of Cortical Receptive Field Profiles," Vision Re- 
search 20, pp. 847-856; and Daugman, J. (1985) "Uncertainty Relation For Resolution In Space, Spatial Frequency, 

25 And Orientation Optimized By Two-Dimensional Visual Cortical Filters," Journal Of The Optical Society Of America, 
Volume 2 (7), pp. 1160-1169. As conjointly optimal filters they offer maximum resolution simultaneously for spatial 
frequency and orientation information, together with 2-D position; they uniquely achieve the lower bound for joint un- 
certainty over these four variables, as dictated by an inescapable uncertainty principle. These properties are particularly 
useful for texture analysis, because of the 2-D spectral specificity as well as positional dependency of texture. 

30 [0036] Two members of the family of 2-D Gabor filters are illustrated in Figure 3, as even-symmetric and odd-sym- 
metric wavelet profiles together with their contour plots. These localized, undulating 2-D functions, defined at many 
different sizes and positions, are multiplied by the raw image pixel data and integrated over their domain of support to 
generate coefficients which describe, extract, and encode image texture information. Applicant has given these the 
name "2-D Gabor filters," because they are a 2-D generalization of a class of elementary functions discussed in one 

35 dimension in 1946 by Dennis Gabor. See Gabor, D. (1946) 'Theory Of Communication," J. Ins. Elec. Eng. , Vol. 93, pp. 
429-457. 

[0037] The 2-D Gabor filters used in the present invention are defined in polar coordinates as follows: 

^ ^ -2*ia>(e-e 0 ) -(r-r Q ) 2 /a 2 -<e-e 0 )V 
40 G(r, 0) = e e e (3) 

where r is radius, 9 is angle in radians, co is angular frequency, and a and (5 are constants. 
[0038] Both the real and imaginary members of the quadrature (even-and odd-symmetric) pair of filters projected 
from the above analytic function are used. Free parameters a and p co-vary in inverse proportion to (o to generate a 
45 multi-scale self-similar family of frequency-selective quadrature filters. They are in quadrature because both orthogonal 
phases are used at each position. They are self-similar because the inverse proportionality of their size and frequency 
parameters renders them ail dilates of each other, sharing a common shape. Their locations, specified by 8 0 and r 0 , 
range across the zones of analysis of the iris. 

[0039] The manner in which an iris code is generated by passing 2-D Gabor filters over the iris, in polar coordinates, 
50 is illustrated in Figures 4A, 4B and 4C. The upper trace (Fig. 4A) shows a 1-D scan around the iris at a particular radius, 

and plots the luminance of the image as a function of the angular coordinate around the iris. (For simplicity, the image 

is represented here as just a 1-D signal rather than a 2-D signal.) The second trace (Fig. 4B) shows the response of 

a Gabor filter with a particular size and symmetry, positioned over each corresponding angular coordinate of the iris. 

It should be noted that because of the bandpass character of Gabor filters, their response to the raw input signal can 
55 be either positive or negative, and is centered around zero. The slow, non-informative shifts in the luminance of the 

raw signal, up and down gradually around the iris resulting from the illumination from below, is removed by the bandpass 

Gabor filters, as is the high frequency noise. 
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[0040] Each bit in an iris code is determined by whether the response of a particular 2-D Gabor filter, having a certain 
size, symmetry, and position over the iris, is positive or negative. This process is denoted in Equations 4, 5, 6 and 7. 
Because it is a "sign bit" that is encoded, this information corresponds to the Most Significant Bit (MSB) of the coefficient 
that results from integrating the product of a 2-D Gabor filter with the input image, as described earlier. The use of both 
5 even and odd quadrature symmetries of 2-D Gabor filters, extracting independent information, is denoted by the sub- 
scripts Re and Im for the bits determined by the corresponding Real and Imaginary parts of the 2-D Gabor filters in 
complex form: 



MSB* e (r, 9) = 1 if Re / / r %ww{i ^e^ r '^e'^ 9 ^I{p t <t>)pdpd<t> > 0 ( 4 ) 
jfi j$ 



« MSB« e (r, 9) = 0 if Re jf j t e -'-('-«) e H'-<)V<>V<'-^'/(p, *)pdpd4> < 0 ( 5 ) 

MSB /m (r,0) = 1 if Im jT jT e -*««(»-*) e -l'-') , /«' ',-(«-•>»/* r(p,<l>)pdpd4 > 0 (6) 

20 



MSB lm (r,6) = 0 if ImJJ^'^^-^e'^ 1 ^*)^ * 0 < 7 > 

25 

[0041] These conditionals (Equations 4-7) determine each of the 2,048 bits in an iris code, across multiple scales of 
analysis (set by parameters a, p and go ) and across all of the sampled positions (set by polar coordinate parameters 
r and 6 ) within the defined zones of analysis of the iris image. It is noteworthy that a very significant amount of data 
30 compression is achieved in such a code because of its decorrelating nature. Whereas the original iris image may 
consist of typically 262,000 bytes (a 512 x 512 array of pixels, each pixel requiring one byte), the significant iris texture 
has been reduced by this multi-scale 2-D Gabor code to a very compact signature comprising only 1/1 ,000th as much 
data (namely 256 bytes). 

[0042] An example of a 256-byte iris code is inscribed at the top of Figure 2, organized as 256 angular columns, 

35 each with 8 bits computed over the concentric zones of analysis. Although there are 2,048 bits in any given code, a 
code possesses fewer than 2,048 independent binary degrees-of-freedom. The chief reason is that there exist sub- 
stantial radial correlations in an iris. For example, a given furrow tends to propagate across a significant radial distance, 
and thus it exerts its influence on several remote parts of the code. A second reason is that correlations are introduced 
by the low-pass component of the bandpass 2-D Gabor filters. Specifically, any signal convolved with a linear filter 

40 acquires a correlation distance equal to the reciprocal of the bandwidth of the filter. 

[0043] The actual number of independent degrees-of-freedom can be estimated by examining the distribution of 
Hamming distances (fraction of disagreeing bits) across a large population of iris codes, comparing each code bit-by- 
bit with every other code computed from a different iris. Since each bit has equal probability of being a 1 or a 0, there 
is a probability p = 0.5 that any pair of bits from different iris codes disagree. If each of the 2,048 bits in a given code 

45 were independent from every other bit, then the distribution of observed Hamming distances should be equivalent to 
a binomial distribution with p = 0.5 and N = 2,048 (in other words, equivalent to tossing a fair coin repeatedly and 
counting the fraction of heads in each round of 2,048 tosses). The actual distribution of observed Hamming distances 
among 2,064 codes from different irises is shown in Figure 5. Its standard e rror is o = 0.038, around a mean of u. = 
0.502. Since the standard deviation of a binomial distribution is given by a = Vpq/N , the observed distribution of Ham- 

50 ming distances is equivalent to a binomial distribution with N = 1 73 bits; an actual fit of such a binomial to the observed 
iris code data, which may be seen in Figure 10, reveals an excellent match. Thus there are approximately 173 inde- 
pendent binary degrees-of-freedom in a 2,048-bit iris code. 

[0044] Using the binomial estimate of N = 173 binary degrees-of-freedom as a measure of the complexity or dimen- 
sionality of a 2,048-bit iris code, we can compute the likelihood of two codes from different irises matching by chance. 
55 Since the 2-D Gabor filters have no positive or negative bias, the a priori odds of any given bit being a 1 or a 0 are 
even, and hence the probability is 0.5 that two corresponding bits in two different iris codes would be the same. Factoring 
in the partial correlations within an iris code but assuming independence between iris codes, the odds that two different 
irises might generate the same code are one in 2 173 , which equals one in 10 52 . 
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[0045] The process of comparing any two iris codes (Figure 1 , Block 26), such as a previously stored one (Block 22) 
and one that is currently computed from a presenting image (Block 24), is very simple because of the universal format 
and fixed length of ail such codes. A similarity metric called a Hamming distance is computed, that measures "distance, 
" or similarity between the two codes. This measure simply adds up the total number of times that two corresponding 

5 bits in the two iris codes disagree. Expressed as a fraction between 0 and 1 , the Hamming distance between any iris 
code and an exact copy of itself would therefore be 0, since all 2,048 corresponding pairs of bits would agree. The 
Hamming distance between any iris code and its complement (in which every bit is just reversed), would be 1. The 
Hamming distance between two random and independent strings of bits would be expected to be 0.5, since any pair 
of corresponding bits has a 50% likelihood of agreeing and a 50% likelihood of disagreeing. Thus, if two iris codes 

10 arise from different eyes, their Hamming distance would be expected to be 0.5; if they arise from the same eye, on 
different occasions, their Hamming distance would be expected to be considerably lower. If both iris codes were com- 
puted from an identical photograph, their Hamming distance should approach zero. 

[0046] Comparisons between iris codes can be made with several different relative shifts along their angular axis, 
in order to compensate for possible tilt of the head of the subject or torsional eye rotation. These relative shifts in the 
15 code comparison process are readily implemented by lateral scrolling of the iris codes relative to each other, as though 
the code portrayed in the upper left corner of Figure 2 were wrapped around into a cylinder, joining the left and right 
margins, and then rotating the cylinder and repeating the comparison process. 

[0047] The computation of Hamming distances between iris codes is made very simple through the use of the ele- 
mentary logical operator XOR (Exclusive-OR). A pair of bits A and B can have exactly four possible combinations: (AB) 
20 = (00), (01), (10), and (11). The XOR operator on two inputs is defined as 1 if one and only one of the inputs equals 1; 
otherwise their XOR is 0. Thus, in the example given for the four possible combinations of the bits A and B, the corre- 
sponding values of their XOR are: (A XOR B) = 0, 1, 1,0. Clearly, XOR can thus be used to detect a disagreement 
between any pair of bits, whatever their values. 

[0048] Tallying up their total number of times that the XOR of two corresponding iris code bits equals 1 , and dividing 
25 by the total number of such comparisons (which is the number of bits in an iris code), is equivalent to measuring the 
Hamming distance between the two codes. Alternatively, this quantity may also be described as the normalized squared 
length, or squared-norm, of the difference vector between the two iris code vectors in a 2,048-dimensional binary 
space. All of these formulations generate the same metric for iris code comparisons, and they lend themselves to 
immediate conversion into a computed probability that two iris codes came from the same iris, and hence from the 
30 same person. 

[0049] The problem of recognizing the signature of a given iris as belonging to a particular individual, or deciding 
that s/he is an imposter, can be formulated within the framework of statistical pattern recognition and decision theory. 
[0050] Yes/No decisions in pattern recognition have four possible outcomes: either a given pattern is or is not a true 
instance of the category in question; and for either of these two cases, the decision made can be correct or incorrect. 

35 These four outcomes are usually termed Hit, Miss, False Alarm, and Correct Rejection. In the present application the 
four possible outcomes are Acceptance of Authentic (AA), Acceptance of Imposter (IA), Rejection of Authentic (AR), 
and Rejection of Imposter (IR). The goal of the decision-making algorithm is to maximize the likeli-hoods of AA and 
IR, while minimizing the likelihoods of IA and AR. The pairwise trade-offs among the probabilities of these four outcomes 
can be manipulated in a way that reflects their associated costs and benefits in a particular application. 

40 [0051] A formulation for decision-under-uncertainty is presented in Figure 6. A given measurement of the Hamming 
distance, or fraction of disagreeing bits between two iris codes, constitutes a point on the abscissa. The measurement 
is regarded as being a random variable describing one of two processes, represented by the two overlapping probability 
distributions. It is unknown a priori which of these two distributions the random variable describes; the goal is to decide 
which. A criterion is chosen, as indicated by the dotted vertical line in Figure 6, and all Hamming distances smaller 

45 than this criterion are judged to belong to the "Authentics" distribution, while all Hamming distances greater than this 
criterion are judged to belong to the "Imposters" distribution. These two distributions, Pa u ( h ) and p imp( H )» 9' ve the 
probability density of a particular measured Hamming distance, H, arising from two comparisons of the same iris (an 
"authentic"), or from two comparisons of different irises (an "imposter"), respectively. 

[0052] The four outcomes AA, I A, AR, and IR have probabilities that are now fully determined by the chosen criterion 
50 and by the statistical parameters of the two underlying distributions. If the decision rule is: 

Accept If Hamming Distance < Criterion 

Reject If Hamming Distance > Criterion 
then the probabilities of the four possible outcomes are equal to the areas under the two probability density functions, 
P Au (H) and P| mp (H), on either side of the chosen criterion, C: 

55 
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P(AA) = J Q °P Au (H)dH 

5 

P{AR) = £ P Au (H)dH 
P(IA) = £ P Im {H)dH 



P(IR) = j l c P Im {H)dH 



20 [0053] These four probabilities are signified by the four shaded areas in Figure 6. 

[0054] It is clear that the four probabilities separate into two pairs which must sum to unity, and two pairs which are 
governed by inequalities: 

„ P(AA) + P(AR) = 1 (12) 



P(/A) + P(/R) = 1 (13) 



P(AA) > P(IA) (14) 



P(IR)>P(AR) (15) 

35 

[0055] It is also clear that the two error rates, P(AR) and P(IA), will be minimized if the two Hamming distance 
distributions, P Au (H) and P| mp (H), have minimal overlap. This can be achieved either by pushing their two means 
farther apart, or by reducing their variances, or both. It should be noted that the two distributions in general will not 
have the same form and variance, as was implied in Figure 6 for simplicity. 
40 [0056] The usefulness, or Identification Power, of a biometric signature method for distinguishing among and recog- 
nizing individuals can be defined in terms of the amount of overlap between these two distributions. Clearly, if there 
were no overlap, it would be possible to make correct decisions 100% of the time. Conversely, the more overlap, the 
higher would be the proportion of errors, regardless of the decision criterion employed. 

[0057] Specifying the decision criterion C in Equations (8)-(11) above provides the ability to elect different decision 
45 strategies most suitable for different applications. For example, in controlling access to bank accounts in Automatic 
Teller Machines, allowing the ARR to be mush higher than 0 may be bad for customer relations, even if it means 
tolerating a higher IAR; after all, the cost of accepting an imposter is at worst the ATM cash withdrawal limit. On the 
other hand, in military or embassy security systems, a far more conservative criterion is demanded, perhaps as stringent 
as only one chance in a million for an imposter to be admitted, even though such a criterion in this case would require 
so tolerating a higher ARR (percentage of authentics rejected for further screening). 

[0058] The manipulation of the decision criterion C, in order to implement different decision strategies, is illustrated 
schematically in Figure 7. A theoretical Authentic Acceptance Rate, or P(AA), is plotted against a theoretical Imposter 
Acceptance Rate, or P(IA), as a locus of points determined by different choices for the decision criterion C as indicated 
in Figure 6. These two figures are only theoretical in order to clarify the nature of the decision problem; as drawn, they 
55 imply far greater uncertainty than actually exists for this biometric recognition system. 

[0059] Equation (14) tells us that the strategy curve shown in Figure 7 will always lie above the diagonal line in this 
probability space. In general, one would like to use decision strategies that generate points as close as possible to the 
upper left corner, since reaching that ideal would mean that all authentics were accepted while all imposters were 



(8) 
(9) 
(10) 

(11) 
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rejected. Clearly, strategies that are excessively conservative or excessively liberal correspond to sliding along the 
curve towards the two diagonal extremes, in which all Subjects, authentics and imposters alike, are either universally 
rejected (lower left) or universally accepted (upper right). Clearly the Identification Power in Figure 7 would thus be 
zero anywhere along the diagonal, and it would equal one in the upper left corner of this space. The overall power of 
5 a detection method, regardless of where one chooses to place the decision criterion along the liberal-conservative 
strategy curve, may therefore be gauged by the length of the line segment joining the diagonal line and the bend in 
the strategy curve, i.e., the length of the "arrow" in the "bow". 

[0060] Having formulated the biometric identification problem within a framework of signal processing and statistical 
decision theory, we can now evaluate the identifiability of persons by their irises. 

10 [0061] The distribution of Hamming distances computed among 2,064 pairs of different irises (labelled "imposters" 
because no pair came from the same person) was seen previously in Figure 5. As expected, the average Hamming 
distance was close to 0.5, since any bit in the 2,048-bit (256-byte) codes for two different irises has equal probability 
of agreeing or disagreeing. The distribution of Hamming distances was tightly clustered around the expected value; 
the actual mean is u. = 0.497 with a = 0.038 standard deviation. 

15 [0062] Figure 8 shows the distribution of Hamming distances computed among 1 ,208 pairs of different images of the 
same iris ("authentics"), obtained at different times. Ideally these Hamming distances should be zero, if the images 
were truly identical; however, differences in angle of gaze, partial eyelid closure, specular reflections from the cornea, 
and relative constriction of the pupil, generate some differences in the structure encoded. Nonetheless, these Hamming 
distances are clearly substantially smaller than those seen in Figure 5 for imposters. The authentics distribution has 

20 a mean of u. = 0.084 with o = 0.0435 standard deviation. 

[0063] Because of the possible variation in the tilt of a person's head on different occasions, together with cyclover- 
gence (torsional rotation) of the eye in its socket, it is necessary to carry out all iris code comparisons over a range of 
different relative orientations. Only the best match obtained from a set of such comparisons is kept as the measure of 
similarity. Because such a "best of n relative orientations" test always selects the lowest Hamming distance, both for 

25 comparisons of authentics and imposters, both of these distributions are shifted towards the left and are made narrower 
than they would be otherwise. Of course, this does not affect the nature of the decision task, which is not based on 
any assumptions about the forms of the two distributions. But it improves overall decision performance, since the 
degree of match between different images of the same eye benefits much more from comparing them at several can- 
didate relative orientations, than does the match between codes for unrelated irises. 

30 [0064] The authentics distribution seen in Figure 8 was obtained with the best of 7 orientations rule. Figure 9 shows 
this distribution in conjunction with the corresponding distribution for imposters, which compares the same set of 2,064 
unrelated iris codes as was used for Figure 5 but now using the same new "best of 7 orientations" rule. These two 
distributions have no empirical overlap. Figure 10 shows the same pair of histograms fitted by theoretical binomial 
distributions, whose parameters are fitted to match the means and variances of the observed pair of empirical distri- 

35 butions. The binomial form is mathematically appropriate, given the nature of the code comparisons process as a 
statistical sequence of Bernoulli (coin-tossing) trials. The problem of iris pattern recognition is thereby transformed 
essentially into a statistical test of independence. The Exclusive-OR test between iris code bits examines the hypothesis 
that the two code sequences could have risen from independent random processes. This test of statistical independ- 
ence is failed almost certainly for two codes arising from the same eye, but it is passed almost certainly for two codes 

40 arising from different eyes. 

[0065] Although the two empirical distributions seen in Figure 9 have no overlap, and indeed there were no observed 
points falling in the range of 0.25 to 0.35 Hamming distance, theoretically there would have been some overlap between 
the two distributions if a large enough database had been studied. The fitted pair of theoretical binomial curves super- 
imposed in Figure 10 provides a way to estimate the error rates if there had been an unlimited number of observations. 

45 These overlapping fitted distributions should be considered with Figure 6 in mind, which introduced the basic framework 
for statistical decision-making. As specified in Equations (8)-(11), the probabilities of personal identity or non-identity, 
and the predicted error rates, can be calculated as the cumulative integrals under these two distributions on either side 
of any chosen Hamming distance decision criterion. 

[0066] The error rates for falsely accepting imposters (IAR) and rejecting authentics (ARR), in the observed popu- 
50 lation of 1,208 paired comparisons, are plotted in Figure 11 as a function of Hamming distance criterion. This graph 
makes explicit the trade-off between Type I and Type II error rates (IAR and ARR), which can be manipulated by 
changes in criterion. Their theoretical cross-over point occurs for a Hamming distance criterion of about 0.321 , at which 
point the probabilities of false acceptance and false rejection are both 1 in 131 ,000. For situations demanding a much 
more conservative decision criterion, such as for example a Hamming distance of 0.26, for which the chances of 
55 imposter acceptance are around one in two billion, Figure 11 shows that the rate of correctly accepting authentics is 
still higher than 99.96% . The dotted curves are the theoretical error rates computed according to Equations (7) and 
(8), using the fitted binomial distributions for P Au (H) and P| mp (H) seen in Figure 10. 

[0067] The binomial framework developed earlier for this statistical pattern recognition task allows us to compute 



11 



EP 0 664 037 B1 



10 



15 



20 



25 



30 



35 



40 



45 



50 



the confidence levels associated with any decision confirming or disconfirming a person's identity on the basis of their 
iris signature. Specifically, when two iris codes are compared, e.g., a previously "enrolled" one and a presenting one, 
and their Hamming distance is measured, we can calculate the odds that a Hamming distance of this size or smaller 
might have arisen by chance from two different irises. Only if this probability is sufficiently small would the person be 
accepted as authentic. 

[0068] The confidence level associated with the decision requires summing over all possible combinations 



of the N independent binary degrees-of-freedom in which there are m < CN bits mismatching (resulting in a Hamming 
Distance HD less than or equal to the fraction C of mismatches), times the probability of each such event, given that 
p is the probability that a pair of bits do not match and q is the probability that they do match: 



[0069] Equation (16) specifies the probability that the Hamming distance HD between codes for different irises might 
happen by chance to be smaller than a given decision criterion C. This lets us determine, for any given decision criterion, 
what would be the odds against chance for a false acceptance. (Stirling's approximation allows estimates of the large 
factorials needed to evaluate Equation 16.) As tabulated in Figure 12, these theoretical odds are roughly one in 2.4 
million for a Hamming distance criterion of 0.30, and the odds rapidly reach "planetary" levels (in the billions) for criteria 
of 0.26 or smaller. Obviously, such improbable errors when using such criteria never occurred in the existing data base. 
Indeed, the imposter acceptance rate observed among the 2,064 iris code comparisons was already at zero for a 
Hamming distance criterion as high as 0.35, at which point the theoretical authentic rejection rate was about one in a 
million (bottom row of Fig. 12). 

[0070] A final gauge of the power of the present biometric signature security system is the computed confidence 
level associated with typical, or average, Hamming distances encountered among different pictures of the same iris. 
As was apparent in the histograms of Figures 8 and 9, the average Hamming distance between two iris codes generated 
from the same iris was 0.084. For this typical case, the confidence level for a decision to accept the individual as 
authentic is truly astronomical. Specifically, the odds that a Hamming distance of this average size or smaller could 
arise from an imposter, according to Equation (16), are one in 10 31 . 

[0071] The analysis performed by this embodiment of the present invention generates an identification code of 256 
bytes of data (2048 bits). This number has been found to yield a reliable identification code, combined with optimized 
processing characteristics. The iris code maps ail different irises into a universal, abstract mathematical code of fixed 
length. This allows code comparisons to be made extremely efficiently, even between relatively "featureless" irises 
(perhaps as the result of low contrast imaging), and those imaged with rich visible texture. 

[0072] The comparison process lends itself directly to simple hardware implementation based on the XOR-gate. Built 
from just three transistors, XOR gates are available in standard semiconductor chips. For example, the standard IC, 
named 74F86, contains four independent XOR gates which can operate at 80 Megahertz, and is available at negligible 
cost. Any comparison between two iris codes, as well as exhaustive searches through large databases of stored iris 
codes, can be implemented extremely rapidly and with inherent parallelism. For example, a circuit board containing a 
32 X 32 array of 74F86 ICs, could exhaustively compare a "presenting" iris code against a population of 160 million 
previously stored iris codes within one second, to establish reliably whether the individual is any one of those persons. 
[0073] Because each bit in the code can be regarded as a binary random variable, the theory of binomial statistics 
has been applied to evaluate the probability of any given fraction of the bits in two different iris codes agreeing just by 
chance. This permits the objective calculation of confidence levels for every decision, using statistical decision theory. 
Insofar as the pattern recognition problem has been converted here into a statistical test of independence on the 
polarities of Gabor coefficients computed over multiple scales of analysis in a dimensionless coordinate system, the 
overall theoretical basis of the present invention can be regarded as a synthesis of unique signal processing methods 
with statistical decision theory. 

[0074] It is understood that the above-described embodiment is merely illustrative of the application. Other embod- 
iments may be readily devised by those skilled in the art which will embody the principles of the invention and fall within 
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the scope thereof. 



Claims 

5 

1. A method for uniquely identifying a particular human being by biometric analysis of the iris of the eye, comprising 
the following steps: 

acquiring (10) an image of an eye of the human to be identified; 
10 isolating and defining (12, 14, 16, 18) the iris of the eye within the image; 

defining (18) a plurality of annular analysis bands within the iris, wherein the radial width of an annular band 
is defined as a fixed fraction of the radial width of the iris; 
analyzing (20) the iris to generate a present iris code; 

comparing (26) said present code with a previously generated reference iris code to generate a measure of 
15 similarity between said present iris code and said reference code; and 

converting (28) said similarity measure into a decision that said iris codes either do or do not arise from the 
same iris. 

2. The method of claim 1 , further comprising the step of: 

20 

calculating (30) a confidence level for the decision. 

3. The method of claim 2, wherein said isolating and defining step further includes the steps of: 

25 defining (12) a circular pupillary boundary between the iris and pupil portions of the image; 

defining (14) another circular boundary between the iris and sclera portions of the image, using arcs that are 
not necessarily concentric with the pupillary boundary; and 

establishing (16) a polar coordinate system on the isolated iris image, the origin of the coordinate system being 
the center of the circular pupillary boundary, 

30 

wherein the radial coordinate is measured as a percentage of the distance between the said circular pupillary 
boundary and said circular boundary between the iris and sclera. 

4. The method of claim 3, wherein: 

35 

said analysis bands exclude certain preselected portions of the iris image likely to be occluded by the eyelids, 
eyelashes, or specular reflection from an illuminator. 

5. The method of claim 4, wherein said analyzing step includes the steps of: 

40 

analyzing the portion of the iris image lying within said annular analysis bands and employing signal processing 
means to generate an iris code for said iris image portion. 

6. The method of claim 5, wherein: 

45 

said signal processing means comprises applying multi-scale, self-similar, two dimensional quadrature band- 
pass filters in polar coordinates to the iris image. 

7. The method of claim 6, wherein: 

50 

said iris code has a fixed number of bits and a universal format for all irises^ 

8. The method of claim 7, wherein said analysis step includes: 

55 applying said bandpass filters to a region of the raw iris image signal to remove luminance bias, to remove 

slow luminance gradients arising from oblique illumination, to remove noise, and to prevent aliasing. 

9. The method of claim 8, wherein: 
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the value of each bit in the iris code is specified as a "1" or "0", by computing over whatever region of the iris 
constitutes the support for the said filters of a given size at a given location, the most-significant-bit of the filter 
outputs according to the following definitions: 



MSB„(r,*) = 1 if Re jT jf e -*~<»-*e-« r -'»''-V< # -*''/(* o)pd pdo > 0 
10 MS8*M) = 0 if Re / # e^ # -«6-« r -*'V«%-<#-«V#»/(^ |#)prf ^ < 0 



MSB /m (r,f) = 1 if lm /^/ < «- 1 ~ l# -^e- ,r -*» 1 ^ , e-« | -«V^/ (Pi ^ )prfpd ^ > 0 
MSIWM) = 0 if lm \J 9 t l ^>-*>e' u -*>*\-l'-" 7 t*f( p%0 ) p dpd0 < 0 

20 

10. The method of claim 9, wherein said comparing step includes the steps of: 

comparing any two iris codes by computing the elementary logical XOR (exclusive-OR logical operation) be- 
tween all their corresponding bits; and 

25 

computing the squared norm of the resulting binary vector; 

wherein this comparison measure is defined as the Hamming distance between the two iris code vectors. 

30 11. The method of claim 1 0, further comprising: 

repeating the comparing step for several different relative shifts of the iris code along its angular axis, to 
compensate for possible tilt of the head of the presenting subject or cyclovergence (i.e., torsional eye rotation). 

35 12. The method of claim 11, wherein said calculating step includes the steps of: 

converting said Hamming distance into a calculated likelihood that the two codes originated from the same 
iris, and hence the same person. 

40 13. The method of claim 12, wherein: 

said calculated likelihood is found by computing the probability that an observed matching fraction of bits in 
the presenting code and reference code could match by chance if the codes were independent, i.e., arising 
from separate irises. 

45 

14. The method of claim 13, wherein: 

the measured Hamming distance is converted into a probability that the two said iris codes are from the same 
eye; 

50 

a preselected criterion is applied to said measured Hamming distance to generate a "yes" or "no" decision; and 
a confidence level for the decision is provided by the calculated probability. 
55 15. The method of claim 14, wherein: 

said circular pupillary boundary is defined by the relationship: 
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w 



where r is the radius of said boundary, Xq and y 0 are center coordinates, and I is image intensity; 
wherein said radius and center coordinates are varied in a preselected pattern. 
16. The method of claim 15, wherein: 



said other circular boundary between the iris and sclera portion of the image is defined by determining the 
distances from said coordinate system origin to the left and right limbus in accordance with the following 
15 relationship: 



20 

where r is the pupillary radius previously defined, 8 is a radial shell distance, l(p,6) is the original image intensity 
in polar coordinates, and <(> equals 0 or n in order to find the limbus in the 3-o'clock or 9-o'clock meridia, 
respectively. 

25 

17. The method of claim 16, wherein: 

said analysis bands include a plurality of annular bands extending completely around said pupil, and a plurality 
of semi-annular bands extending from polar angles of approximately 45 degrees to 1 35 degrees and from 225 
30 degrees to 31 5 degrees around said pupil. 

18. The method of claim 17, wherein: 

said analysis bands are spaced at equal fractional radial distances from an inner point located at a preselected 
35 fractional distance from said pupillary boundary to an outer point located at preselected fractional distance 

from said limbus, said spacing being angularly weighted to account for said difference between fractional 
distances from said coordinate center to said right and left limbi. 

19. The method of claim 14, wherein: 

40 

said bandpass filters are 2-D Gabor filters. 

20. The method of claim 19, wherein: 

45 said 2-D Gabor filters are defined in polar coordinates as follows: 

50 

where r is radius, 0 in angular distance in radians, co is frequency, and a and p are constants. 

21. The method of claim 20, wherein: 

55 said iris code is comprised of 2048 bits. 

22. The method of claim 3, wherein: 
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said circular pupillary boundary is defined by the relationship: 



where r is the radius of said boundary, Xq and y Q are center coordinates, and I is image intensity; 
wherein said radius and center coordinates are systematically varied by iterative gradient ascent to find the 
maximum in said defined relationship. 

23. The method of claim 22, wherein: 

said other circular boundary between the iris and sclera portion of the image is defined by determining the 
distances from said coordinate system origin to the left and right limbus in accordance with the following 
relationship: 



where r 0 is the pupillary radius previously defined, 5 is a radial shell distance, l(p,6) is the original image 
intensity in polar coordinates, and § equals 0 or n in order to find the limbus in the 3-o'clock or 9-o'clock meridia, 
respectively. 

24. Apparatus for uniquely identifying a particular human being by biometric analysis of the iris of the eye, comprising: 

imaging equipment for acquiring a digital image of an eye of the human to be identified; and 
a processor, wherein said processor is programmed to isolate and define the iris of the eye within the image, 
define a plurality of annular analysis bands within the iris, with the radial width of an annular band being defined 
as a fixed fraction of the radial width of the iris, analyse the iris to generate a present iris code, compare said 
present code with a previously generated reference iris code to generate a measure of similarity between said 
present iris code and said reference code and convert said similarity measure into a decision that said iris 
codes either do or do not arise from the same iris. 



Pate ntansp ruche 

1 . Verfahren zum eindeutigen Identifizieren einer bestimmten menschlichen Person durch biometrische Analyse der 
Iris des Auges, mit den folgenden Schritten: 

Beschaffen (10) eines Bildes eines Auges der zu identifizierenden menschlichen Person, 
Isolieren und Definieren (12, 14, 16, 18) der Iris des Auges innerhalb des Bildes, 

Definieren (18) einer Mehrzahl ringformiger Analysebander innerhalb der Iris, wobei die radiate Breite eines 
ringformigen Bandes als fester Bruchteil der radialen Weite der Iris festgelegt wird, 
Analysieren (20) der Iris, urn einen aktuell vorliegenden Iriscode zu erzeugen, 

Vergleichen (26) des vorliegenden Codes mit einem fruher erzeugten Referenzcode der Iris, urn ein Mafc der 
Ahnlichkeit zwischen dem vorliegenden Iriscode und dem Referenzcode zu erzeugen, und 
Umwandeln (28) des Mafies der Ahnlichkeit in eine Entscheidung, dali die Iriscodes entweder von derselben 
Iris stammen oder nicht. 

2. Verfahren nach Anspruch 1, welches weiterhin den Schritt aufweist: 

Berechnen (30) eines Vertrauensniveaus (Konfidenz) fur die Entscheidung. 

3. Verfahren nach Anspruch 2, wobei der Schritt des Isolierens und Definierens weiterhin die Schritte aufweist: 
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Definieren einer kreisformigen Pupillengrenze zwischen den Iris- und den Pupillenabschnitten des Bildes, 
Definieren (14) einer weiteren kreisformigen Grenze zwischen den Iris- und den Hornhaut-Abschnitten des 
Bildes unter Verwendung von Bogen bzw. Kreisbogen, die nicht notwendigerweise konzentrisch zu der Pupil- 
lengrenze liegen, und 

5 Bereitstellen (16) eines Systems von Polarkoordinaten fur das isolierte Irisbild, wobei der Ursprung des Ko- 

ordinatensystems das Zentrum der kreisformigen Pupillengrenze ist, wobei die radiale Koordinate gemessen 
wird als Prozentwert des Abstandes zwischen der kreisformigen Pupillengrenze und der kreisformigen Grenze 
zwischen der Iris und der Lederhaut. 

10 4. Verfahren nach Anspruch 3, wobei die Analysebander bestimmte vorgewahlte Abschnitte des Irisbildes ausschlie- 
Ren, die wahrscheinlich durch die Augenlider, die Wimpern oder spiegelnde Reflexion von einer Beleuchtung ver- 
deckt bzw. unsichtbar gemacht werden. 



15 
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5. Verfahren nach Anspruch 4, wobei der Analyseschritt die Schritte aufweist: 

Analysieren des Teiles des Irisbildes, welches innerhalb der ringformigen Analysebander liegt und Verwenden 
von Signalverarbeitungseinrichtungen, um einen Iriscode fur diesen Irisbildbereich zu erzeugen. 

6. Verfahren nach Anspruch 5, wobei 

die Signalverarbeitungseinrichtung das Anwenden von mehrskaligen, selbst-ahnlichen, zweidimensionalen 
Bandpalifiltern mit Phasenverschiebung in Polarkoordinaten fur das Irisbild aufweist. 

7. Verfahren nach Anspruch 6, wobei der Iriscode eine feste Anzahl von Bits und ein universelles Format fur alle 
25 Irisbilder hat. 

8. Verfahren nach Anspruch 7, wobei der Analyseschritt aufweist: 

Anwenden der Band pa (ifi Iter auf einen Bereich des rohen Irisbildsignales, um einen Luminanzfehler zu be- 
30 seitigen, um geringformige Luminanzgradienten zu beseitigen, die von einer schragen Beleuchtung herruhren, 

um Rauschen zu entfernen und um einen Treppeneffekt (Sagezahnverzerrung) zu vermeiden. 

9. Verfahren nach Anspruch 8, wobei: 

35 der Wert jedes Bits in dem Iriscode als "1" Oder "0" bestimmt wird, indem uber irgendeinen Bereich der Iris, 

welcher die Grundlage fur die Filter auf einer gegebenen Grofte und an einer gegebenen Stelle bildet, das 
hdchstwertige Bit (MSB) der Filterausgange gemafi den folgenden Definitionen berechnet wird: 

40 MSB*c(r y 6) = 1, falls Re £ ^^^e^^^^^Iip, ftp dpd<f>>0 

MSB**(r, 9) = 0, falls Re £ ^-"•^« e -<^) 2/ « , e^> l, ^/(^ l ftp dpdj><0 
M55im(r, 9) = 1. falls Im £ ^ 7 ^^^ ia9 e^^ i ^I(p 9 ftp dpd<f»0 

MSB im(r, 9) = 0, falls Im y-"«w e *r-P?'*' e -^M*I(p 9 ftp dpd<f><0 

55 10. Verfahren nach Anspruch 9, wobei der Vergleichsschritt die Schritte aufweist: 

Vergleichen irgendwelcher zwei Iriscodes durch Berechnen des elementaren logischen XOR (logische Aus- 
schliefilich-ODER-Operation) zwischen all ihren entsprechenden Bits, und 
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Berechnen der quadratischen Norm des resultierenden binaren Vektors, 
wobei dieses Vergleichsmafi als der Hamming-Abstand zwischen den beiden Iriscodevektoren definiert wird. 

11. Verfahren nach Anspruch 10, welches weiterhin aufweist: 

Wiederholen des Schrittes fur mehrere verschiedene relative Verschiebungen des Iriscodes entlang ihrer Win- 
kelachse, um eine mogliche Verkippung bzw. Verdrehung des Kopfes des dargebotenen Objektes oder eine 
Zyklovergenz (das heiflt eine Augenverdrehung mit Torsion) auszugleichen. 

12. Verfahren nach Anspruch 11, wobei der Berechnungsschritt die Schritte aufweist: 

Umwandeln des Hamming-Abstandes in eine berechnete Wahrscheinlichkeit, dafi die beiden Codes aus der- 
selben Iris und damit von derselben Person stammen. 



13. Verfahren nach Anspruch 12, wobei 

die berechnete Wahrscheinlichkeit erhalten wird durch Berechnen der Wahrscheinlichkeit, daft ein beobach- 
teter, passender Bruchteil von Bits in dem vorliegenden Code und dem Bezugscode zufallig miteinander uber- 
20 einstimmen konnten, falls die Codes unabhangig voneinander waren, das heifit von verschiedenen Irissen 

herruhren wurden. 

14. Verfahren nach Anspruch 13, wobei: 

25 der gemessene Hamming-Abstand in eine Wahrscheinlichkeit umgewandelt wird, dafi die beiden Iriscodes 

von demselben Auge stammen, 

ein vorgewahltes Kriterium auf den gemessenen Hamming-Abstand angewendet wird, um eine "ja"- oder 
"nein"-Entscheidung zu erzeugen, und 

durch die berechnete Wahrscheinlichkeit ein Konfidenzniveau fur die Entscheidung geliefert wird. 



15. Verfahren nach Anspruch 14, wobei: 

die kreisformige Pupillengrenze definiert wird durch die Beziehung: 

d 



max(r,x 0 .>-o)||— j r ^ ^ /(*, y)ds\\ 



wobei r der Radius der Grenze, Xq und y 0 Koordinaten des Zentrums und I die Intensitat des Bildes sind, 
40 wobei die Koordinate des Radius und des Zentrums nach einem vorgewahlten Muster variiert werden. 

16. Verfahren nach Anspruch 15, wobei: 

die andere kreisformige Grenze zwischen dem Iris- und dem Hornhautabschnitt des Bildes definiert wird durch 
45 Bestimmen der Abstande von dem Ursprung des Koordinatensystems zu dem linken und rechten Bindehaut- 

ring gemafi der folgenden Beziehung: 



d rr+S 2 



wobei r der zuvor definierte Pupillenradius ist, 8 der Abstand bzw. das Mafi einer Radialschale ist, l(p,9) die 
55 ursprungliche Bildintensitat in Polarkoordinaten ist und <(> gleich 0 oder n ist, um die Bindehautgrenze in den Me- 

ridianbereichen der 3 Uhr- bzw. 9 Uhr-Position zu finden. 

17. Verfahren nach Anspruch 16, wobei: 
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die Analysebander eine Mehrzahl ringformiger Bander umfassen, die sich vollstandig um die Pupille herum 
erstrecken, sowie eine Mehrzahl von halbkreisformigen Bandern, die sich von Polarwinkeln von 45° bis 135° 
und von 225° bis 315° um die Pupille herum erstrecken. 

5 18. Verfahren nach Anspruch 17, wobei: 

die Analysebander in gleichen Bruchteilen des radialen Abstandes von einem inneren Punkt aus, der unter 
einem vorgewahlten Bruchteilabstand von der Pupillengrenze angeordnet ist, bis zu einem aufieren Punkt, 
der unter einem vorgewahlten Bruchteilabstand von dem Bindehautring angeordnet ist, voneinander beab- 
10 standet sind, wobei der Abstand in Winkelrichtung gewichtet wird, um den Unterschied zwischen den Bruch- 

teilabstanden von dem Koordinatenzentrum zu den rechten und linken Bindehautringen zu berucksichtigen. 

19. Verfahren nach Anspruch 14, wobei die Bandpafifilter 2-D-Gaborfilter sind. 

15 20. Verfahren nach Anspruch 19, wobei die 2-D-Gaborfilter folgendermalien in Polarkoordinaten definiert sind: 

-2«/cD<e-e 0 ) -(r-r 0 ) 2 /a 2 -<e-e 0 )V 
G(r,8) = e e e 

20 wobei r der Radius, 9 der Winkelabstand in Radian, co die Frequenz und a und p Konstanten sind. 

21. Verfahren nach Anspruch 20, wobei der Iriscode aus 2048 Bits besteht. 

22. Verfahren nach Anspruch 3, wobei die kreisformige Pupillengrenze durch die Beziehung definiert wird: 



25 



30 



40 



45 



wobei r der Radius der Grenze, Xq und y 0 die Koordinaten des Zentrums und I die Bildintensitat ist, 
wobei der Radius und die Zentrumskoordinaten systematisch durch einen iterativen Gradientenanstieg va- 
riiert werden, um das Maximum in der definierten Beziehung zu finden. 

35 23. Verfahren nach Anspruch 22, wobei: 

die andere kreisformige Grenze zwischen der Iris und dem Lederhautabschnitt des Bildes definiert wird durch 
Bestimmen der Abstande vom Ursprung des Koordinatensystems zu der linken und rechten Bindehautgrenze 
gemafi der folgenden Beziehung: 



d rr+S 2 f *wr/8 

max Hi,5r 0 ,, 0 , 0 | y r W*/ 8 l ^ 6)p dp de 

wobei r 0 der zuvor definierte Pupillenradius ist, 8 ein Radialschalenabstand ist, l(p,8) die ursprungliche Bildin- 
tensitat in Polarkoordinaten ist, und § gleich 0 Oder ic ist, um die Bindehautgrenze an dem der Position 3 Uhr bzw. 
9. Uhr entsprechenden Meridian zu finden. 

so 24. Vorrichtung fur das eindeutige Identifizieren einer bestimmten menschlichen Person durch biometrische Analyse 
der Iris des Auges, mit: 

einer Abbildungsausrustung fur das Erfassen eines digitalen Bildes eines Auges der Person, die es zu iden- 
tifizieren gilt, und 

55 einem Prozessor, wobei der Prozessor so programmiert ist, daft er die Iris des Auges innerhalb des Bildes 

isoliert und festlegt, eine Mehrzahl von ringformigen Analysebandern innerhalb der Iris definiert, wobei die 
radiale Weite eines ringformigen Bandes als fester Bruchteil der radialen Weite der Iris definiert ist, weiterhin 
die Iris analysiert, um einen Code der vorliegenden Iris zu erzeugen, den voriiegenden Code mit einem fruher 
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erzeugten Bezugscode der Iris vergleicht, um ein Mad der Ahnlichkeit zwischen dem vorliegenden Iriscode 
und dem Bezugscode zu erzeugen und um das Ahnlichkeitsmafi in eine Entscheidung umzuwandeln, daft die 
Iriscodes entweder von derselben Iris herriihren oder nicht. 



Revendications 

1. Procede pour identifier de facon unique un etre humain particulier par analyse biometrique de Tiris de Toeil, com- 
prenant les etapes suivantes : 

10 

I'acquisition (10) d'une image d'un oeil de Tetre humain a identifier; 

I'isolation et la definition (12, 14, 16, 18) de Tin's de I'oeil a I'interieur de Timage ; 

la definition (18) d'une pluralite de bandes d'analyse annutaires a I'interieur de I'iris, la largeur radiale d'une 
bande annulaire etant definie comme etant une fraction fixe de la largeur radiale de I'iris ; 
is Panalyse (20) de I'iris pour generer un code d'iris present ; 

la comparaison (26) dudit code present a un code d'iris de reference genere precedemment pour generer une 
mesure de ressemblance entre ledit code d'iris present et ledit code de reference ; et 
la conversion (28) de ladite mesure de ressemblance en une decision selon laquelle lesdits codes d'iris pro- 
viennent ou ne proviennent pas du meme iris. 

20 

2. Procede selon la revendication 1 , comprenant en outre I'etape : 

de calcul (30) d'un niveau de confiance pour la decision. 

25 3. Procede selon la revendication 2, dans lequel ladite etape d'isolation et de definition comprend en outre les etapes 
consistant a : 

definir (12) une limite pupillaire circulate entre les parties formant i'iris et la pupille de Timage ; 
definir (14) une autre limite circulaire entre les parties formant I'iris et la sclerotique de Timage, en utilisant des 
arcs qui ne sont pas necessairement concentriques avec la limite pupillaire ; et 

etablir (16) un systeme de coordonnees polaires sur Timage d'iris isolee, Torigine du systeme de coordonnees 
etant le centre de la limite pupillaire circulaire, les coordonnees radiales etant mesurees comme etant un 
pourcentage de la distance entre ladite limite pupillaire circulaire et ladite limite circulaire entre I'iris et la scle- 
rotique. 

Procede selon la revendication 3, dans lequel : 

lesdites bandes d'analyse excluent certaines parties preselection nees de Timage d'iris susceptibles d'etre 
occultees par les paupieres, les cils ou une reflexion speculate provenant d'un illuminateur. 

40 

5. Procede selon (a revendication 4, dans lequel ladite etape d'analyse comprend les etapes consistant a : 

analyser la partie de Timage d'iris qui se situe a I'interieur desdites bandes d'analyse annulaires, et a utiliser 
des moyens de traitement du signal pour generer un code d'iris pour ladite partie d'image d'iris. 

45 

6. Procede selon la revendication 5, dans lequel : 

ledit moyen de traitement du signal comprend Tapplication a Timage d'iris de filtres passe-bande bidimension- 
nels en quadrature, auto-similaires et a echelles multiples en coordonnees polaires. 

50 

7. Procede selon la revendication 6, dans lequel : 

ledit code d'iris a un nombre fixe de bits et un format universel pour tous les iris. 

55 8. Procede selon la revendication 7, dans lequel ladite etape d'analyse consiste a : 

appliquer lesdits filtres passe-bande a une region du signal d'image d'iris brut pour eliminer le biais de lumi- 
nance, pour eliminer les faibles gradients de luminance resultant d'un eclairement oblique, pour eliminer le 
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bruit et pour empecher un repliement spectral. 

9. Procede selon la revendication 8, dans lequel : 

la valeur de chaque bit dans le code d'iris est specifiee par un "1" ou un "0", en effectuant un calcul dans toute 
region de i'iris qui constitue un support pour lesdits filtres d'une taille donnee a un emplacement donne, le bit 
de poids le plus fort des sorties du filtre etant tel que defini ci-apres : 

MSB*,(r t *) = 1 if Re ^/^-^-•»e-''"'"''* , e-« | -^^/(p i *)/Hfp^ > 0 
MSB*,M) = 0 if RtfJ 4 c -^^^l c -U-.) J /o^^^/^7(p^)pdp^ < 0 

MSB, m (r.*) = 1 if lm/J^" 2 ~ t '-' , e-^^°%- (, - d,V ^V(p^)p^ > 0 

MSB, m M) = 0 \r\m IJ 9 t ^-"t- u -*''*\-V-* )%i **I{p.<p)pdpdo < 0 

25 10. Procede selon la revendication 9, dans lequel ladite etape de comparaison comprend les etapes consistant a : 

comparer deux codes d'iris quelconques en calculant le OU-X logique elementaire (operation logique OU 
exclusif) entre tous leurs bits correspondants ; et 
calculer la norme elevee au carre du vecteur binaire resultant ; 
30 cette mesure de comparaison etant definie comme etant la distance de Hamming entre les deux vecteurs de 

codes d'iris. 

11. Procede selon la revendication 10, consistant en outre : 

35 a repeter I'etape de comparaison pour plusieurs decalages relatifs differents du code d'iris le long de son axe 

angulaire, afin de compenser I'eventuelle inclinaison de la tete du sujet se presentant ou une cyclovergence 
(c'est-a-dire une rotation oculaire de torsion). 
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12. Procede selon la revendication 11, dans lequel ladite etape de calcul comporte les etapes consistant a : 

convertir ladite distance de Hamming en une vraisemblance calculee selon laquelle tes deux codes prove- 
naient du meme iris et par consequent, de la meme personne. 

13. Procede selon la revendication 12, dans lequel : 

ladite vraisemblance calculee est trouvee par calcul de la probabilite selon laquelle une fraction concordante 
observee des bits dans le code de presentation et dans le code de reference peut concorder par hasard si 
les codes sont independants, c'est-a-dire s'ils proviennent d'iris distincts. 

50 14. Procede selon la revendication 13, dans lequel : 

la distance de Hamming mesuree est convertie en une probabilite selon laquelle lesdits deux codes d'iris 
proviennent du meme oeil ; 

un critere preselectionne est applique a ladite distance de Hamming mesuree pour generer une decision de 
55 type "oui" ou "non" ; et 

un niveau de confiance pour ta decision est fourni par la probabilite calculee. 

15. Procede selon la revendication 14, dans lequel : 
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ladite limite pupillaire circulate est definie par la relation suivante : 



ou r est le rayon de ladite limite, Xq et y 0 sont des coordonnees centrales, et I est I'intensite de I'image ; 

dans lequel ledit rayon et lesdites coordonnees centrales sont amenees a varier selon un motif preselec- 
tionne. 

16. Procede selon la revendication 15, dans lequel : 

ladite autre limite circulate entre les parties formant I'iris et la sclerotique de I'image est definie en determinant 
les distances de ladite origine du systeme de coordonnees vers les limbes gauche et droit conformement a 
la relation suivante : 



oil r est le rayon pupillaire precedemment defini, 8 est une distance d'enveloppe radiale, l(p, 9) est I'intensite 
d'image originale en coordonnees polaires, et $ est egal a 0 ou n afin de trouver respectivement le limbe suivant 
les meridiens de 3 heures ou de 9 heures. 

17. Procede selon la revendication 16, dans lequel : 

lesdites bandes d'analyse incluent une pluralite de bandes annulaires s'etendant entierement autour de ladite 
pupille, et une pluralite de bandes semi-annulaires s'etendant depuis des angles polaires d'environ 45 degres 
a 135 degres et de 225 degres a 315 degres autour de ladite pupille. 

18. Procede selon la revendication 17, dans lequel : 

lesdites bandes d'analyse sont espacees a des distances radiales fractionnaires egales d'un point interieur 
situe a une distance fractionnaire preselectionnee de ladite limite pupillaire a un point exterieur situe a une 
distance fractionnaire preselectionnee dudit limbe, ledit espacement etant pondere angulairement pour tenir 
compte de ladite difference entre des distances fractionnaires dudit centre de coordonnees auxdits limbes 
droit et gauche. 

19. Procede selon la revendication 14, dans lequel : 

lesdits filtres passe-bande sont des filtres de Gabor bidimensionnels. 

20. Procede selon la revendication 19, dans lequel : 

lesdits filtre de Gabor bidimensionnels sont definis en coordonnees polaires de la facon suivante : 



ou r est le rayon, 9 est la distance angulaire en radians, to est la frequence, et a et (5 sont des constantes. 
21. Procede selon la revendication 20, dans lequel : 
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ledit code d'iris est constitue de 2048 bits. 



22. Procede selon la revendication 3, dans lequel : 



ladite limite pupillaire circulaire est definie par la relation suivante : 



ou r est le rayon de ladite limite, Xq et y 0 sont des coordonnees centrales, et I est I'intensite de I'image ; 
ledit rayon et lesdites coordonnees du centre etant amenes a varier systematiquement par une remontee 
iterative suivant le gradient, pourtrouver le maximum de ladite relation definie. 

23. Procede selon la revendication 22, dans lequel : 

ladite autre limite circulaire entre les parties formant flris et la sclerotique de I'image est definie en determinant 
les distances de ladite origine du systeme de coordonnees aux limbes gauche et droit conformement a la 
relation suivante : 



ou r 0 est le rayon pupillaire precedemment defini, 8 est une distance d'enveloppe radiale, l(p,0) est I'intensite 
d'image originale en coordonnees polaires, et § est egal aOoun afin de trouver le limbe selon les meridiens 
respectifs de 3 heures ou de 9 heures. 

24. Appareil destine a identifier de fagon unique un etre humain particulier par analyse biometrique de I'iris de I'oeil, 
comprenant : 

un equipement de formation d'image destine a I'acquisition d'une image numerique d'un oeil de I'etre humain 
a identifier ; et 

un processeur, ledit processeur etant programme pour isoler et definir I'iris de I'oeil a I'interieur de I'image, 
definir une pluralite de bandes d'analyse annulaires a I'interieur de I'iris, la largeur radiale d'une bande annulaire 
etant definie comme etant une fraction fixe de la largeur radiale de I'iris, analyser I'iris afin de generer un code 
d'iris present, comparer ledit code present a un code d'iris de reference genere precedemment afin de generer 
une mesure de ressemblance entre ledit code d'iris present et ledit code de reference, et convertir ladite 
mesure de ressemblance en une decision selon laquelle lesdits codes d'iris proviennent ou ne proviennent 
pas du meme iris. 
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| ACQUIRE IMAGEj ^'O 



DEFINE PUPILLARY IRIS BOUNDARY] ^ 12 



DEFINE LIMBIC IRIS BOUNDARY.L^' 
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ESTABLISH COORDINATE SYSTEM] ^" 



DEFINE ANALYSIS BANDS] ^' 
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ANALYZE IMAGE DATA-[ ^~ 
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STORE REFERENCE CODE 



PRESENT IDENTIFICATION CODEj ^ 24 



COMPARE CODES TO COMPUTE HAMMING DISTANCES] 
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IDENTIFY OR REJECT SUBJECT^ 
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CALCULATE DECISION CONFIDENCE LEVEL 



-30 



FIG. I 
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FIG. 2 
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FIG. 3 A 



FIG. 3B 
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FIG. 4A 
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